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We utilize the fluid-gravity duality to investigate large order behavior of hydrodynamic gradient 
expansion of the dynamics of a gauge theory plasma system. This corresponds to the inclusion of 
dissipative terms and transport coefficients of very high order. Using the dual gravity description, 
we calculate numerically the form of the stress tensor for a boost-invariant flow in a hydrodynamic 
expansion up to terms with 240 derivatives. We observe a factorial growth of gradient contributions 
at large orders, which indicates a zero radius of convergence of the hydrodynamic series. Further- 
more, we identify the leading singularity in the Borel transform of the hydrodynamic energy density 
with the lowest non-hydrodynamic excitation corresponding to a 'non-hydrodynamic' quasinormal 
mode on the gravity side. 



Introduction. Hydrodynamics is an effective theory of 
fluids describing long-wavelength evolution of conserved 
quantities, among which are the energy (local tempera- 
ture) and momentum (local velocity). In the relativistic 
setting of uncharged hydrodynamics, which we consider 
here, one starts with the perfect fluid stress tensor and 
systematically corrects it by including all available gra- 
dient terms (dissipative corrections) with their relevance 
decreasing with the number of derivatives they carry. 
Each gradient term enters with an a priori distinct func- 
tion of temperature - transport coefficient, e.g. the first 
order symmetric traceless contribution to the stress ten- 
sor is associated with the shear viscosity. 

Transport coefficients are, in principle, derivable from 
the microscopic theory underlying the given fluid, being 
it kinetic theory for fluids formed out of weakly coupled 
quantum fields or the gauge-gravity duality for certain 
strongly coupled liquids. The exact values of transport 
coefficients are of course important for the understanding 
of the dynamics of the particular fluid. A famous exam- 
ple is the shear viscosity of quark-gluon plasma created 
at RHIC and LHC experiments, which needs to be suffi- 
ciently small in order to account for the observed spectra 
of particles (see [T] for a recent review). 

Given the ubiquity and universality of the hydrody- 
namic description of many physical systems, it is in- 
teresting to investigate the character of the hydrody- 
namic gradient expansion and its interrelation with non- 
hydrodynamic collective excitations/degrees of freedom. 

On a formal level, it is important to understand 
the large-order behavior of any approximate scheme in 
physics based on the existence of a small parameter such 
as hydrodynamic gradient expansion. Given the com- 
plexity of typical perturbative calculations at large or- 
ders, the question we are asking is challenging enough to 
deserve attention and interest. In fact, no methods were 
available so far to proceed effectively beyond the first few 
terms in the hydrodynamic gradient expansion. 



The issues of high order hydrodynamics are also in- 
teresting from a phenomenological point of view, in the 
context of ongoing relativistic heavy ion collision pro- 
grams [T. Two examples considered recently in the lit- 
erature are refining the criterium for the applicability of 
hydrodynamics by including terms with more than two 
derivatives [21 13] and trying to resum hydrodynamic ex- 
pansion in order to account for situations in which indi- 
vidual gradient terms are not small enough to truncate 
the series at lowest orders [H [S] . 

In this letter we provide strong evidence that the ra- 
dius of convergence of the hydrodynamic gradient expan- 
sion is zero. Furthermore, using the standard, in the 
context of divergent series, technique of Borel transform, 
we identify the leading singularity in Borel-transformed 
hydrodynamic stress tensor with the lowest lying non- 
hydrodynamic degree of freedom. As a byproduct, we un- 
cover an intriguing structural similarity of the phenom- 
ena arising in perturbative expansions in coupling con- 
stant in quantum systems with temporal evolution gov- 
erned by hydrodynamic and non-hydrodynamic modes. 

As a way of generating hydrodynamic stress ten- 
sor at high orders of gradient expansion, we utilize 
the AdS/CFT correspondence [B] and the fluid-gravity 
duality [7 . We obtain gradient-expanded stress ten- 
sor for a particular solution of the equations of (all- 
order) relativistic hydrodynamics, generalizing the so- 
called Bjorken flow [HI [9] - a toy model of expanding 
quark-gluon plasma. The underlying microscopic the- 
ory is strongly coupled N — A super Yang-Mills at large 
number of colors (A^c) and at strong coupling. The low- 
est lying non-hydrodynamic degree of freedom is then the 
metric quasinormal mode sharing the symmetries of the 
problem and having the lowest frequency. 
Bjorken flow and holography. Bjorken flow is a par- 
ticular solution of the equations of relativistic hydrody- 
namics and describes matter expanding in one dimension 
with the additional assumption of boost-invariance in the 
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longitudinal direction. This symmetry can be made man- 
ifest upon passing to curvilinear proper time r - rapidity 
y coordinates related to the lab frame time x'^ and posi- 
tion along the expansion axis via 

x'^ = T coshy and x^ — t sinhy. (1) 

In the case of (3-|-l)-dimensional conformal field theory 
plasma, the most general stress tensor obeying the sym- 
metries of the problem in coordinates (t, y, a;^, a;^) reads 

n = diag(-e,pi,pT,PT)^., (2) 

where the energy density e is a function of proper time 
only and the longitudinal and transverse px pressures 
are fully expressed in terms of the energy density [9j 

PL = —f — Te' and pT = £+-^Te'. (3) 

Note that in the proper time - rapidity coordinates ([T]) 
there is no momentum flow in the stress tensor ^ and 
so the flow velocity is trivial and takes the form u — dr- 
Hydrodynamic constituent relations lead then to gradient 
expanded energy density of the form 

' = (A$7^ ['^ + ^3/(A-)'^' + ^4/(Ar)^/3 + ...), 

(4) 

where £2 = 1 and A is the only free parameter setting an 
overall energy scale. The prefactor was chosen to match 
Af — 4: super Yang-Mills at large- A^c and strong coupling. 
In the following, without a loss of generality we set A = 

Large-r expansion of the energy density in powers of 
r^^/'', as in Q, is equivalent to the hydrodynamic gra- 
dient expansion and arises from expressing gradients of 
velocity (V^Uj^ ~ t^^) m units of the effective tempera- 
ture (T ~ e^/^ ~ T~^^^). The value of the coefficient ei 
is related to the shear viscosity 77, whereas £2 is a sum 
of two transport coefficients: relaxation time rn and the 
so-called Ai [10]. Higher order contributions to the en- 
ergy density are expected to be linear combinations of so 
far unidentified transport coefficients. Note also that the 
expansion Q is sensitive to both linear and nonlinear 
gradient terms. 

As explained in [11] [12] , higher order contributions to 
the energy density Q can be obtained by solving Ein- 
stein's equations with negative cosmological constant for 
the metric ansatz of the form 

ds^ = 2dTdr-AdT^ + j:^e~^^dy^+^^e^{dxl+dxl), (5) 

where the warp factors A, S and B are functions of r 
and T constructed in the gradient expansion as required 
by the fluid-gravity duality. At leading order, the warp- 
factors are that of a locally boosted black brane and this 
solution gets systematically corrected in t"^/^ expansion, 



as is the case with the energy density in the dual field 
theory Q. 

The background expanded in t"^/'^ around locally 
boosted black brane is slowly evolving and captures 
only hydrodynamic degrees of freedom. One can in ad- 
dition consider the incorporation of non-hydrodynamic 
(fast evolving) degrees of freedom by linearizing Ein- 
stein's equations on top of the hydrodynamic solution, 
i.e. B — Bfiyiij-o + SB, and similarly for A and S, and 
looking for SB corresponding to (at very large time) ex- 
ponentially decaying contribution to the stress tensor 
depending only on r. For the static background anal- 
ogous calculation would lead to the spectrum of non- 
hydrodynamic quasinormal modes carrying zero momen- 
tum, which is known to be the same as the spectrum 
of zero momentum quasinormal modes for the massless 
scalar field |T5] . 

In the leading order of the gradient expansion, the re- 
sulting modes, on the gravity side, indeed essentially re- 
duce to the scalar quasinormal modes but obtain an ad- 
ditional factor of | and are damped exponentially in ts 
[14) . Upon including viscous correction, the modes ob- 
tain a further nontrivial power-like pre-exponential factor 

Se - t"'"" exp {-i - LUqnrrl T^^^)- (6) 

Explicit gravity calculation for the lowest mode yield 
w^„„ = 3.1195-2.7467, a,„,„ = -1.5422-f 0.5199 i. (7) 

The frequency Wg„„i agrees with the frequency of the 
lowest non-hydrodynamic scalar quasinormal mode and 
was calculated before in [14], whereas the prediction of 
c^qnm IS a new result specific to the dissipative modifi- 
cations of the expanding black hole geometry. In the 
following we will be able to reproduce numerically ([t]) 
just from the large order behavior of the hydrodynamic 
series. 

Large order behavior of hydrodynamic energy 
density. Numerical implementation of the methods out- 
lined in [TTJ [T2] allow for efficient calculation of hydrody- 
namic series given by Q, up to a very large order, since 
one is effectively solving a set of linear ODE's (coming 
from Einstein's equations) at each order. Using spectral 
methods we iteratively solved these equations in the large 
time expansion reconstructing the energy density up to 
the order 240, i.e. up to the term £242 in To the best 
of our knowledge this is the first approach allowing to 
access information about the large order behavior of the 
hydrodynamic series in any physical system or model. 

As a way of monitoring the accuracy of our procedures 
we compared normalized values of evaluated Einstein's 
equations at each order of the t"^/^ expansion to the 
ratio of coefficients of gradient-expanded energy density 
to gradient expanded warp factors. This ensures that 
our results for the energy density are reliable. We also 
verified that we reproduced the known analytic results 
for the energy density at low (< 3) orders. 
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FIG. 1. Behavior of the coefficients of hydrodynamic series 
for the energy density as a function of the order. At large 
enough order the coefficients start exhibiting factorial growth. 
The radius of convergence of the Borel transformed series is 
estimated to be 6.37, in rough agreement with ([9|: 2/3 x 
6.37 = 4.25. 



As anticipated in the introduction, the coefficients of 
the gradient expanded energy density Q that we ob- 
tained, indicate that the hydrodynamic gradient expan- 
sion, as seen in the example of boost-invariant flow, has 
zero radius of convergence. This is clearly visible on 
Fig. [T] which shows the factorial growth with order of 
the coefficients of the hydrodynamic series for the energy 
density Q. 

Singularities of the Borel transform and quasinor- 
mal modes. A standard way of dealing with divergent 
series, including perturbative quantum field theories [T5] . 
is performing a Borel transform of the original power se- 
ries (here a series in m = t~3), e„ = e„/ n\, and subse- 
quently the Borel resummation 
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(8) 



where the integral is taken over the real positive axis. 
A key issue in performing Borel resummation is the iden- 
tification of the singularities of the Borel transform and 
its analytical continuation into some neighborhood of the 
positive real axis. The singularities of the Borel trans- 
form are interesting for their own sake as typically they 
have a definite physical interpretation (like instanton 
contributions in quantum field theories). We will find 
that this will also be the case in our context. 

Since the Borel transform has typically only a finite 
radius of convergence, we use the standard technique of 
Pade approximants to provide an analytic continuation. 
Such an approach was used with success in the context of 
resumming perturbative series in quantum field theories 
in |16j and, in principle, can be refined by including the 
information about the behavior of the series at large val- 
ues of the parameter (typically strong coupling behavior 
for perturbative expansion in quantum field theories). 

In our case this regime corresponds to small times 
(high order dissipative terms). Although it was estab- 
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FIG. 2. Real and imaginary parts of poles (o of the symmet- 
ric Pade approximant of the Borel transform of the energy 
density (4|. From the plot we removed numerically spurious 
poles. The pole closest to the origin governs the convergence 
radius of the Borel transformed series and gives rise to the 
lowest quasinormal frequency. The poles from the encircled 
region (magenta) lead to the power-like pre-exponential factor 
in (|6|. 



lished that, when expanded around t = 0, the energy 
density contains only even powers of proper time |17j . 
this information turns out to be hard to implement in an 
unambiguous way and in the following we adopted the 
simplest analytic continuation. 

Fig. [5] shows the position of poles of the symmetric 
Pade approximant of the Borel transformed energy den- 
sity. Let us note that some care must be taken as the 
Pade approximant exhibits apparent poles on the nega- 
tive real axis which however cancel almost perfectly (up 
to 10"^"° accuracy) with zeroes of the numerator. As our 
knowledge about the series is limited to a finite number 
of terms (first 241), we should only trust the structure 
in the close vicinity of the origin. Note however that 
the lack of poles on the positive real axis seems to in- 
dicate Borel-summability, pointing towards the possible 
existence of Borel-resummed all-order hydrodynamics. 

The poles approximate some complicated structure of 
branch-cuts. The pole nearest the origin, from which a 
major branch-cut starts, sets the radius of convergence 
of the Borel transform of hydrodynamic series. Its nu- 
merical value (multiplied by the factor 3/2i) reads 

WBo,,, = 3.1193 - 2.747H (|cjBo,e;| = 4.1565) (9) 

and is consistent with the coefficient obtained from fitting 
factorial-like behavior of the original series. 

Deforming the contour from the real positive axis to 
the one encircling some of the zeros of the denomina- 
tor of the Pade approximated Borel transform, through 
residues, leads to an ambiguity in resummation. It is 
tempting to speculate that such an ambiguity can be un- 
derstood as the appearance of new, non-hydrodynamic, 
degrees of freedom in the system. Indeed, perform- 
ing contour integration around the cut originating from 
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(3/2i) X uJBorei leads to the following contribution to the 
energy density at large proper time 



Se 



where 



OiBorel 



exp {-i 3/2 UJBorel T^^^), 



-1.5426 + 0.5192i. 



(10) 



(11) 



This might lead to refined criterium for the applicability 
of hydrodynamics used in [2l [3] and goes very much in 
the spirit of [U 



The contribution (10) together with ([9| and (11) 
matches, up to numerical accuracy, the gravitational pre- 
diction for the behavior of the lowest metric quasinormal 
mode given by (|6]) and ([?]), but is derived entirely within 
the hydrodynamic gradient expansion. 
Summary and conclusions. We used the fluid-gravity 
duality to numerically obtain the form of hydrodynamic 
stress tensor of the Bjorken flow up to 240*'* order in gra- 
dients. This corresponds to the inclusion of dissipative 
terms of a very high order. We discovered that the cor- 
rections grow factorially with order. This provides strong 
indication that hydrodynamic expansion is an example of 
asymptotic series, i.e. it has a zero radius of convergence. 

Upon a simple analytic continuation of Borel trans- 
formed series we found that it is the frequency of the 
lowest non-hydrodynamic quasinormal mode which sets 
the radius of convergence of the Borel transform. Subse- 
quently, from the hydrodynamic expression alone we ex- 
tracted the modification of this contribution due to the 
shear viscosity and matched with the gravity calculation. 
It is fascinating to observe that non-hydrodynamic col- 
lective excitations, which are not describable within the 
framework of hydrodynamics nevertheless leave their im- 
print in the behavior of high order transport coefficients. 
This is reminiscent of the relation of instanton contri- 
butions with the divergence of the perturbative series in 
quantum field theory. 

Note that our observation about the asymptotic char- 
acter of the hydrodynamic series is made upon evaluating 
the stress tensor of a particular type of hydrodynamic 
flow. The high degree of symmetry of that flow - boost- 
invariance - means that at each order in the large proper 
time expansion, we are indeed studying new dissipative 
terms in the hydrodynamic stress tensor with new trans- 
port coefficients. The dependence on initial conditions 
for this hydrodynamic flow reduces just to a trivial over- 
all rescaling. 

We speculate that the most likely source of the facto- 
rial growth of the coefficients of the hydrodynamic series 
for the energy density is the fast growth of the number of 
gradient terms contributing to the hydrodynamic stress 
tensor at each order, in analogy with the factorial growth 
of the number of Feynman graphs at large order of per- 
turbative calculations. 

Regarding future directions, a thought-provoking phe- 
nomenological spin-off of our letter is the possibility of 
resumming hydrodynamic description and extending the 
hydrodynamic stress tensor past the regime in which sub- 
sequent low order terms give comparable contributions. 
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